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HIGHLIGHTS 


►  We  analyzed  experimentally  water  transport  and  flooding  in  DMFC. 

►  We  proposed  and  validated  a  model  of  water  transport  and  flooding  in  DMFC. 

►  Water  flow  through  cathode  diffusion  layer  is  due  to  vapor  diffusion  and  liquid  permeation. 

►  Flooding  implies  both  superficial  and  bulk  pores  obstruction  of  diffusion  layer. 

►  We  developed  correlations  to  simulate  flooding  effects. 
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Water  transport  control  is  the  major  issue  of  direct  methanol  fuel  cell.  High  water  flow  rates  through  the  fuel 
cell  imply  extra  water  feeding  at  the  anode  and  flooding  at  the  cathode.  In  the  present  work,  water  transport 
and  flooding  in  the  direct  methanol  fuel  cell  are  investigated  through  both  experimental  and  modeling 
analyses  and  an  interpretation  of  such  phenomena  is  proposed.  The  model  is  validated  on  the  experimental 
data  of  two  different  fuel  cells  in  an  extensive  range  of  operating  conditions.  The  analysis  elucidates  that  water 
transport  through  the  cathode  diffusion  layer  is  determined  by  vapor  diffusion,  slightly  affected  by  current 
density,  and  by  liquid  water  permeation  proportional  to  current  density,  that  occurs  when  liquid  pressure  in 
the  electrode  exceeds  a  threshold  value.  To  simulate  the  effects  of  cathode  diffusion  layer  flooding  two 
mechanisms  must  be  considered  simultaneously:  superficial  pore  obstruction,  proportional  to  liquid  water 
concentration  in  cathode  channel,  and  bulk  pore  obstruction,  proportional  to  liquid  water  permeation.  The 
modeling  analysis  proposes  the  correlations  to  reproduce  the  effects  of  cathode  flooding  and  permits  to 
discuss  the  onset  and  magnitude  of  such  phenomenon  and  the  influence  of  micro-porous  layer. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  receives  considerable 
attention  as  a  leading  candidate  power  source  for  portable  and 
automotive  applications.  The  widely  use  of  the  DMFC  technology  is 
still  hindered  by  some  technological  issues,  among  which  the  water 
management  is  one  of  the  most  important  [1—4],  DMFCs  are  fed 
with  a  liquid  mixture  composed  by  water  and  methanol  and  the 
water  as  well  as  the  methanol  flow  through  the  membrane  from 
the  anode  to  the  cathode  [5,6].  Water  crossover  through  the 
membrane  may  cause  two  problems  named  anodic  water 
consumption  [7-9]  and  cathode  flooding  [10-17],  The  increased 
water  consumption,  caused  by  water  crossover  through  the 
membrane,  implies  the  necessity  of  additional  water  feeding  and 
storage,  reducing  the  compactness  of  the  DMFC  systems  [10,11], 
Secondly,  the  high  rate  of  water  crossover  increases  the  cathode 
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flooding,  especially  during  dynamic  operation,  therefore  a  suitable 
management  of  water  transport  and  flooding  is  considerably 
important  in  order  to  improve  DMFC  performance  and  lifetime. 

The  present  work  aims  to  increase  the  understanding  of  water 
transport  and  flooding  in  DMFC,  proposing  a  model  able  to  describe 
and  reproduce  experimental  data  and  an  exhaustive  interpretation 
of  such  phenomena. 

During  the  last  years,  water  crossover  in  DMFCs  has  been 
extensively  studied  [10—24],  It  has  been  found  that  the  cathode  gas 
diffusion  layer  (GDL),  composed  by  a  carbon  macro-porous  layer 
and  usually  by  an  additional  micro-porous  layer  (MPL)  coated  on  it, 
plays  an  important  role  in  both  water  crossover  through  the 
membrane  and  water  removal  to  the  cathode.  The  effect  of  cathode 
GDL  properties  on  water  transport  in  DMFCs  has  been  widely 
investigated  [10-17,22-25],  Several  works  studied  also  how  to 
minimize  water  transport,  especially  in  passive  DMFC,  modifying 
the  GDLs  characteristics  in  anode  and  cathode  [17,26-28],  The 
characterization  of  flooding  onset  and  its  effect  on  the  performance 
has  not  been  significantly  examined,  despite  this  phenomenon 
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could  occur  also  in  optimized  DMFC  during  dynamic  operation  or 
after  GDL  degradation  [29],  Modeling  analyses  dealt  with  water 
transport  in  DMFC  [30-40],  Schultz  and  Sundmacher  [33]  pre¬ 
sented  a  one-dimensional  dynamic  model,  validated  also  on  water 
transport  experimental  measurements,  neglecting  the  two-phase 
behavior  in  both  diffusion  layers  and  flooding  effects.  Ko  et  al. 
[34]  developed  a  1-D  model  where  the  two-phase  behavior  was 
considered,  but  water  diffusion  in  gas  phase  through  the  cathode 
GDL  was  neglected  and  the  model  was  validated  on  a  single  oper¬ 
ating  condition.  Yang  and  Zhao  [35,36]  presented  a  complete  two- 
dimensional  and  two-phase  DMFC  model,  validated  only  on  the 
experimental  performance  and  a  limited  set  of  operating  condi¬ 
tions  was  considered.  In  previous  works  [37-39],  Casalegno  et  al. 
presented  a  1-D  model  validated  with  respect  to  overall  polariza¬ 
tion,  anode  polarization,  methanol  cross-over,  showing  a  good 
agreement  with  the  experimental  data,  but  water  transport  was  not 
analyzed  in  detail.  Pasaogullari  and  Wang  [40]  developed  a  two- 
phase  numerical  model  to  analyze  the  multi-layer  GDL  in 
PEMFCs.  They  particularly  investigated  the  effects  of  porosity, 
thickness  and  hydrophobicity  of  an  MPL  on  the  two-phase  trans¬ 
port  [40],  The  results  showed  that  the  addition  of  the  MPL  between 
catalyst  layer  and  GDL  could  enhance  liquid  water  removal  and 
reduce  the  liquid  saturation  in  the  catalyst  layer. 

In  the  literature,  the  mechanisms  that  regulate  water  transport 
in  DMFC  are  not  fully  consolidated  and  limited  effort  has  been 
dedicated  to  analyze  flooding  onset  and  its  consequences  on  DMFC 
operation.  A  detailed  comprehension  of  such  phenomena  is 
required  to  further  optimize  components  and  operation  strategies. 

This  work  aims  to  propose  a  combined  experimental  and 
modeling  analysis  in  order  to  increase  the  understanding  of  the 
mechanisms  that  regulate  water  transport  and  flooding  in  DMFC 
through  the  following  steps: 

•  a  model  of  water  transport  and  cathode  flooding  is  proposed 
analyzing  experimental  data  obtained  with  different  operating 
conditions  and  cathode  GDL  features; 


the  previously  developed  DMFC  model  [37-39]  is  integrated, 
in  order  to  validate  the  proposed  interpretation  in  an  extensive 
range  of  operating  conditions,  using  three  different  sets  of 
measure:  polarization  curves,  methanol  crossover  and  water 
flow-rate  at  cathode  outlet; 

correlations  able  to  reproduce  the  effect  of  cathode  GDL 
flooding  are  proposed  as  a  tool  for  further  studies  and 
components  optimization. 


2.  Experimental  setup 

The  experimental  analyses  of  DMFC  performance  and  methanol 
crossover  are  carried  out  implementing,  in  the  equipment  previ¬ 
ously  developed  [41],  the  measure  of  water  concentration  at 
cathode  outlet. 

Cathode  side  was  modified  introducing,  Fig.  1 : 

•  a  heat  exchanger  to  warm  up  the  cathode  exhaust  to  evaporate 
eventual  liquid  water; 

•  a  thermo-hygrometer  for  humidity  (uncertainty  2%)  and 
temperature  (uncertainty  0.2  °C)  measurements,  located  in 
a  thermo-stated  housing  to  avoid  condensation. 

The  above  modifications  permit  to  measure  the  water  content  in 
the  cathode  exhaust  Xf,  0  and  considering  the  mass  balance 
equations  (l)-(8)  reported  in  Ref.  [41],  the  water  flow  at  the 
cathode  outlet  is  given  by: 

Nout°C  =  XH20'Nouttot  (1) 

These  new  components  have  negligible  pressure  drop  and  do  not 
increase  considerably  the  transient  time  necessary  to  reach  steady 
state  in  the  vessels. 

The  experimental  analysis  is  carried  out  on  two  of  the  three 
MEAs  already  characterized  in  Ref.  [41  ] :  MEA  GG,  GDL  without  MPL 


Fig.  1.  Experimental  setup  scheme. 
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Table  1 

Investigated  operating  conditions. 
MEA  T  [K]  Xch,oh  [wt%] 
GM  333  325 

GM  333  3.25 

GM  333  3.25 

GM  333  3.25 

GM  333  6.5 

GM  333  6.5 

GM  333  6.5 

GM  333  6.5 

GM  353  3.25 

GM  353  3.25 

GM  353  6.5 

GM  353  6.5 

GG  333  3.25 

GG  333  3.25 

GG  333  6.5 

GG  333  6.5 

GG  333  6.5 

GG  353  6.5 


maj  [g  min  ]] 


1.14 

1.14 

0.62 

1.14 

0.62 

1.14 


Pin  [kPa]  pout  [kPa] 
112  108 

151  147 

122  115 

155  148 

112  108 

151  147 

122  115 

155  148 

112  108 

122  115 

112  108 

122  115 

112  108 

122  115 

112  108 

122  115 

155  148 

122  115 


on  both  anode  and  cathode;  MEA  GM,  GDL  without  MPL  at  anode 
and  GDL  with  MPL  at  cathode. 

A  set  of  the  measurement  already  examined  in  Ref.  [41  ]  have 
been  reproduced  by  means  of  the  improved  experimental  setup. 
Note  that  in  all  the  investigated  operating  conditions  (Table  l1)  the 
anode  flow  rate  is  set  at  1  g  min-1  and  the  anode  mean  pressure  at 
105  kPa.  The  results  of  performance  and  methanol  crossover 
reported  in  Ref.  [41],  as  well  as  measurement  uncertainty  and 
reproducibility  evaluation,  remain  valid.  The  uncertainty  associated 
to  the  measurement  of  water  flow  is  evaluated,  according  to  Ref. 
[42],  equal  to  7%  in  the  whole  investigated  range. 


3.  Experimental  analyses 

3.1.  Influence  of  operating  conditions 

In  Fig.  2  the  water  flow  at  cathode  outlet,  in  function  of  current 
density,  is  reported  for  different  operating  conditions.  A  general 
behavior  can  be  distinguished:  at  low  current  density  the  water 
flow  is  approximately  constant,  a  sort  of  plateau  is  evident,  then  it 
increases,  in  some  cases  suddenly,  with  an  approximately  linear 
trend.  This  behavior  suggests  the  presence  of  two  different  mech¬ 
anisms  regulating  water  transport  to  cathode  channel.  A  possible 
interpretation  could  be  based  on  the  two  main  mechanisms  of 
water  transport  through  GDL  reported  in  literature:  diffusion  in 
vapor  phase  and  permeation  in  liquid  phase.  The  plateau  could  be 
caused  by  water  vapor  diffusion  through  the  cathode  GDL,  deter¬ 
mined  by  an  approximately  constant  water  concentration  along  the 
cathode  electrode.  Small  current  density  variations  do  not  affect 
significantly  water  concentration  profile.  Instead  at  higher  current 
density  the  linear  trend  could  be  due  to  the  onset  of  liquid  water 
permeation,  forced  by  enhanced  electro-osmotic  drag  and  water 
production  at  cathode,  both  directly  proportional  to  current 
density. 

This  preliminary  interpretation  is  discussed  analyzing  the 
experimental  results  in  different  operating  conditions. 

The  proposed  plateau  origin  is  coherent  with  the  following 
considerations  regarding  Fig.  2a: 

•  an  increase  of  airflow  rate  induces  a  decrease  of  water 
concentration  in  cathode  channel,  as  a  consequence  the  water 


1  In  the  followings  unpressurized  cathode  is  named  as  1  bar,  instead  pressurized 
cathode  is  named  as  1.5  bar. 


concentration  3.25%;  b)  cathode  pressure  1  bar. 

concentration  difference  across  the  GDL  and  the  relative 
diffusive  flow  grow;  therefore  the  water  flow  at  cathode  outlet, 
i.e.  the  plateau  value,  increases; 

•  an  increase  of  temperature  determines  a  general  increment  of 
both  GDL  diffusivity  and  water  concentration  at  cathode  elec¬ 
trode,  related  to  saturation  concentration;  thus  the  water  flow 
at  cathode  outlet  increases; 

•  at  higher  cathode  mean  pressure,  the  total  gas  mixture 
concentration  increases,  therefore  at  similar  water  concentra¬ 
tion,  water  fraction  0  lowers;  consequently  the  plateau 
value  decreases,  according  to  equation  (1). 

Regarding  the  linear  behavior  at  higher  current  density,  the 
curves  at  lower  temperature  assume  the  same  linear  trend, 
coherently  with  the  proposed  dependence  on  electro-osmotic  drag 
and  water  production,  not  affected  considerably  by  the  airflow  or 
pressure  variation. 

A  further  confirmation  of  the  proposed  interpretation  is  obtained 
comparing  the  experimental  results  at  higher  methanol  concen¬ 
tration,  Fig.  2b.  Doubling  methanol  concentration,  the  plateau  values 
do  not  change  significantly,  while  the  performances  are  consider¬ 
ably  different,  Fig.  3,  as  well  as  the  methanol  crossover  [41  ].  There¬ 
fore,  the  plateau  is  confirmed  to  be  influenced  exclusively  by  cathode 
flow  conditions.  Moreover  a  linear  trend  with  a  similar  slope  is 
present,  as  expected  in  these  operating  conditions,  Fig.  2b. 

An  additional  evaluation  of  the  proposed  interpretation  is 
reported  in  the  modeling  section. 
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A  consequence  of  the  considerable  water  transport  through  the 
cathode  GDL  is  the  pore  flooding  that  hinders  oxygen  transport  and 
DMFC  performance.  Liquid  water  permeation  through  the  GDL  is 
expected  to  cause  a  significant  pores  obstruction.  A  comparison  of 
polarization  curves  in  different  operating  conditions  is  reported  in 
Fig.  3.  At  low  methanol  concentration  an  increase  of  airflow  rate 
determines  negligible  variations  of  DMFC  performance.  In  contrast, 
water  flow  at  cathode  outlet  varies  significantly  and  considerable 
water  permeation  is  expected  to  occur  in  both  cases  at  high  current 
density.  This  behavior  can  be  explained  considering  that  the 
performance  at  high  current  density  is  strongly  limited  by  the 
anode  feeding,  worsen  by  low  methanol  concentration,  as  dis¬ 
cussed  previously  in  Ref.  [41  ].  In  such  conditions,  the  effect  of  an 
oxygen  transport  hindered  by  flooding  is  undetectable. 

The  behavior  is  different  at  higher  methanol  concentration:  the 
performance  is  not  anymore  limited  by  anode  feeding  but  a  nega¬ 
tive  effect  of  airflow  reduction  is  evident.  The  water  transport 
increases  at  high  current  density  when  the  airflow  decreases, 
Fig.  2b.  Both  the  lower  airflow  rate  and  the  increased  water  flow 
imply  a  reduction  of  oxygen  concentration  at  cathode  electrode 
leading  to  a  decrease  in  the  performance.  The  performance 
decrease  could  be  due  to  such  oxygen  concentration  reduction  or  to 
a  further  pore  flooding  effect,  but  it  cannot  be  distinguished  from 
these  experimental  results.  A  further  interpretation  of  this  aspect  is 
discussed  in  the  next  section. 


3.2.  Cathode  MPL  influence 

In  order  to  clarify  the  effect  of  water  transport  on  DMFC  perfor¬ 
mance,  the  results  regarding  the  MEA  GG  (without  cathode  MPL) 
and  the  MEA  GM  (with  cathode  MPL)  are  discussed  and  compared.  In 
Fig.  4a,  the  water  flow  at  the  cathode  outlet  for  the  two  MEAs  is 
reported.  An  enhanced  water  transport  is  evident  in  the  MEA 
without  MPL:  the  mass  transfer  resistance  at  cathode  side  decreases 
eliminating  the  additional  micro-porous  layer,  coherently  with  the 
results  reported  in  Ref.  [43],  The  plateau  rises  approximately  by 
a  factor  1.5  in  both  the  two  temperature  conditions.  Moreover  it  is 
confirmed  that  the  enhanced  methanol  crossover  is  due  to  the 
increment  of  water  transport  that  causes  an  increase  of  the  average 
methanol  concentration  at  anode  side,  as  suggested  in  Ref.  [41  ]. 

A  general  improvement  of  performance  is  evident  with  cathode 
without  MPL,  Fig.  4b,  accordingly  to  the  enhanced  mass  transport 
through  the  GDL  that  increases  oxygen  concentration  at  cathode 


electrode.  A  clear  exception  is  evident:  the  performance  of  MEA  GG 
drops  and  becomes  lower  than  MEA  GM  one  at  current  density 
higher  than  0.35  A  cm-2.  This  behavior  has  been  already  found  in 
literature  [44—46]  and  it  is  probably  caused  by  a  considerable  GDL 
flooding  phenomenon. 

In  order  to  confirm  flooding  presence  and  quantify  its  effect 
a  specific  modeling  analysis  is  reported  in  the  followings. 

4.  Modeling 

In  this  work  the  model  already  presented  in  Refs.  [37,38]  is 
integrated  with  equations,  that  describe  water  transport 
phenomena  and  flooding  effects.  The  original  model  includes  the 
complete  mass  balance,  considering  two-phase  anode  flow,  but 
single-phase  cathode  flow  and  consequently  cathode  flooding  was 
not  considered.  Moreover  the  water  crossover  is  calculated  as  the 
sum  of  electro-osmosis  and  liquid  phase  diffusion,  neglecting  liquid 
water  concentration  at  cathode  side. 

In  the  new  model  the  implementation  of  water  transport 
phenomena,  including  two-phase  cathode  flow  and  flooding,  is 
introduced  to  validate  the  proposed  interpretations  and  elucidate 
flooding  effects  on  performance. 

4.J.  Water  transport  description 

According  to  the  proposed  interpretation,  the  water  fluxes 
through  a  DMFC  are  those  illustrated  in  Fig.  5.  The  water  availability 
at  anode  side  is  considered  to  be  very  high  and  the  anode  GDL  is 
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Fig.  5.  Water  fluxes  through  the  DMFC. 

assumed  to  be  flooded  with  fully  liquid  pathways  [17],  therefore  the 
water  transport  through  the  MEA  is  generally  regulated  by  water 
transport  through  membrane  and  cathode  diffusion  layer,  in 
agreement  with  the  literature  [44], 


widely  accepted  in  literature  [40],  Generally,  under  typical  oper¬ 
ating  conditions  p{ijL0  is  greater  than  p^jL0,  leading  to  a  permeation 
flux  from  cathode  to  anode.  Both  anode  and  cathode  total  pressures 
are  assumed  to  vary  linearly  with  channel  length  x: 

Pin  -  (pfn  -  Pom)  *  =  Pa’c  (7) 

The  flux  of  water  due  to  liquid  diffusion  can  be  expressed  as: 


where  and  represent,  respectively,  liquid  water 

concentrations  at  the  membrane  interfaces.  Generally,  the  magni¬ 
tude  of  both  C”£L  and  C™  q,l  depends  on  the  properties  of  the  MEA 
and  the  operating  conditions.  Moreover  in  literature  magnitude 
and  direction  of  water  diffusion  flux  are  controversial  [50,51],  In 
this  work  C™2o’L  is  assumed  equal  to  the  liquid  concentration  cor¬ 
responding  to  a  fully  hydrated  membrane,  with  water  content 
about  30  [52];  Cj^L0  is  already  defined  in  equation  (6). 

Equation  (2)  substitutes  equation  (27)  of  the  original  model  [37], 
The  variable  P^0  is  determined  considering  the  water  transport 
through  gas  diffusion  layer. 


4.1.1.  Membrane 

The  water  cross-over  through  the  membrane  is  due  to  three 
transport  mechanisms  [14,47,48],  electro-osmotic  drag  by  proton 
transport,  diffusion  by  water  concentration  gradient  and  convection 
by  hydraulic  pressure  gradient,  as  reported  in  the  following  equation: 


The  water  flux  due  to  electro-osmotic  drag  can  be  determined  from: 


N* 


i2o 

i,drag 


(3) 


Permeation  flux  depends  on  the  liquid  pressure  difference  across 
the  membrane  and  can  be  expressed  as: 


_ Xjn  Pi  hO 

^H20'MH20-lm 


(Ph2Lo-Ph2Lo) 


(4) 


Since  the  flow  velocity  through  the  liquid  pathways  of  anode  GDL 
can  be  neglected,  pJjL0  is  assumed  equal  to  the  total  pressure  in  the 
anode  flow  channel  pa.  At  cathode  side  material  hydrophobicity  and 
water  transport  could  determine  a  considerable  capillary  pressure, 
thus  pji£0  can  be  significantly  greater  than  gas  pressure.  The  capillary 
pressure  is  proportional  to  water  saturation  in  the  electrode  [49]  and 
consequently  to  liquid  concentration  at  membrane  interface  [27]: 


Ape  =  Ph2L0  -  Ph'To^C.o1  (5) 

This  dependence  is  uncertain  and  strongly  affected  by  material 
characteristics,  in  fact  a  general  correlation  is  not  available  in  liter¬ 
ature.2  In  the  present  work  a  first  order  dependence  is  assumed: 


ch2oL  =  w2-(p£0-P^o)  (6) 

where  W2  is  calibrated  over  experimental  data. 

Due  to  negligible  gas  pressure  gradient  across  the  entire  porous 
media,  p^j®0  can  be  assumed  equal  to  total  cathode  pressure  pc,  as 


4.1.2.  Cathode  diffusion  layer 

The  water  flow  through  cathode  GDL,  N^°,  is  equal  to  the  sum 
of  water  cross-over  through  the  membrane,  N^Soss.  and  water 
sources  in  the  cathode  catalyst  layer  due  to  methanol  cross-over 
oxidation,  ox,  and  to  oxygen  reduction,  according  to: 

*§2  =  WCH%,ox  +  <2,?d  +  <2c°ross  O) 


(10) 


As  reported  in  Ref.  [53],  the  water  flow  through  GDL  is  due  to 
different  contributions:  diffusion  and  permeation  in  gas  phase  and 
permeation  in  liquid  phase.  The  assumption  of  a  constant  gas- 
phase  pressure  across  the  entire  porous  media  [40]  allows 
neglecting  gas  permeation  through  the  GDL.  Therefore  water 
transport  through  cathode  GDL  occurs  mainly  by  vapor  diffusion, 
Wcdl  diff ■  and  liquid  permeation,  perm,  according  to  the 
proposed  interpretation. 

The  two-phase  behavior  at  cathode  side  is  simplified  as  a  plug 
flow,  where  plugs  and  bubbles  have  the  same  local  velocity. 
Equation  (14)  of  the  original  model  [37]  is  replaced  by  the 
following: 


(11) 


C^2o  is  the  time-averaged  total  water  concentration  in  the  cathode 
channel,  defined  as: 


(Ch20  =  C^g0,  if Ch2q  <  C^q 

Ch2o  =  e6 '  q, ?o  +  (1  ■ -  ■ O  -|g|,  if  CcH2o>C^0 

where  f  is  the  cathode  volumetric  gas  fraction. 

According  to  Ref.  [43],  the  overall  mass  transport  coefficient  in 
gas-phase  can  be  determined  from: 


2  Note  that  the  correlations  reported  in  literature  [27]  are  generally  between  the  pG  _  pG  ,  pC,  _  1  ,  (  IgDL  ,  JmPlA 

saturation  and  the  dissolved  water  concentration  in  the  ionomer  of  the  catalyst  H2°  H20,conv  H20,diff  hconv  \Dgdl  '  DmPL/ 

layer,  that  is  different  from  membrane  one,  C™2oL. 


(13) 
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Fig.  6.  a)  Superficial  pores  obstruction,  b)  Bulk  pores  obstruction. 


where  hconv  is  the  convective  mass  transport  coefficient  in  the 
channel.  The  resulting  water  flux  can  thus  be  expressed  as: 


(Ct,c 

y-H2Q  n20 ) 

4n 


(14) 


where  Cf£0,  the  water  vapor  concentration  in  the  cathode  elec¬ 
trode,  depends  on  the  liquid  concentration  in  the  membrane 
interface  C^  q1";  the  data  reported  in  Ref.  [54]  are  interpolated  with 
a  third  order  polynomial. 

The  interpretation  of  water  permeation  onset  at  high  current 
density  is  reported  in  the  experimental  analysis  discussion  and  it  is 
implemented  in  the  model  as  follows.  Water  diffusion  through 
cathode  GDL  led  to  an  increase  of  water  concentration  along  the 
cathode  channel:  this  transport  mechanism  becomes  less  and  less 
considerable  along  cathode  channel  when  its  driving  force 
diminishes. 

This  implies  an  increase  of  the  liquid  pressure  in  the  cathode 
electrode,  p[jL0:  when  it  becomes  greater  than  a  threshold  value, 
Plim,  liquid  permeation  through  cathode  GDL  occurs.  The  values  of 
pf£0  and  Ngdl  perm  are  calculated  by  solving  the  following  water 
balance  at  electrode-GDL  interface: 


N 


f  AP* 


—  IV"2 
r  ,vm. 


[perm  +  NCH3OH,ox 
'GDL.diff  _l“  '''GDL, perm 

considering  the  following  two  cases 

N! 


fN 


02,rd 


=  NHl 


ifPH,o>Pi™ 


'GDL, perm  —  u 

N™  °erm=  (Plim_PH2o) 

^Sdiff  =D^2°- (Ch2oLW2  -p,im) 

=nJ12~ 


v  '"GDL, perm 

+WH2° 

+,',ch3oh, 


m.drag-1- 


•+<2 


-NHl 


oxygen  transport  caused  by  GDL  and  electrode  pore  obstruction 
with  liquid  water  [34,36,49],  The  following  model  does  not  pretend 
to  be  a  detailed  and  exhaustive  description  of  such  complex 
phenomenon,  but  a  tool  to  increase  its  understanding  and  for 
further  studies.  In  this  work,  flooding  effects  are  modeled  consid¬ 
ering  two  contributions: 

•  superficial  pores  obstruction,  Fig.  6a:  if  water  concentration  in 
cathode  channel  exceeds  the  saturation  value,  drops  of  water 
condense  on  GDL  surface.3  This  effect  is  considered  as 
a  reduction  of  GDL  diffusivity  proportional  to  liquid  water 
concentration  [43],  where  correlations  for  similar  GDLs  are 
proposed.  A  suitable  one  is  adapted  as  follows: 

DGDL,02  =  DGDL,eff,02  “  ’  (CH20  “  CH2o)  07) 

where  B i  is  derived  from  the  model  calibration,  Section  4.2; 

•  bulk  pores  obstruction,  Fig.  6b:  if  liquid  water  permeation 
through  cathode  GDL  occurs,  the  flow  obstructs  GDL  pores, 
establishing  liquid  pathways  [55],  Both  effects  imply  a  reduc¬ 
tion  of  GDL  diffusivity  proportional  to  liquid  water  permeation 
through  cathode  GDL.4  To  take  into  account  this  effect 
a  correlation  is  proposed: 

^GDL,02  =  CGDL,eff,02  _  B2  •  (WGDL,perm)  (18) 

where  B2  and  B3  are  obtained  by  model  calibration. 


4.3.  Model  validation 

The  proposed  DMFC  model  is  constituted  by  7  differential 
equations,  29  algebraic  equations,  36  variables  and  two  correlations 
to  take  into  account  flooding  effects,  equations  (17)  and  (18).  This 
DAE  system  is  solved  in  MATLAB®  environment,  applying  the 
appropriate  initial  conditions,  regarding  inlet  flows,  potential 
difference  and  fuel  cell  temperature  and  pressures.  The  values  of 
the  parameters  utilized  for  the  calculation  are  reported  in  Table  2. 

The  calibration  procedure  consists  in  the  minimization  of  the 
residuals  between  model  estimation  and  experimental  results  and 
permits  to  determine  the  values  of  the  most  uncertain  physic- 
chemical  parameters.  The  considered  experimental  data  set  for 
calibration  concern  the  MEA  GM  and  is  composed  of  132 
measurement  points,  coming  from  12  polarization  curves,  Table  1. 


4.2.  Flooding  description 

It  is  widely  known  that  flooding  reduces  fuel  cell  performance, 
while  the  effective  mechanisms  are  not  fully  understood.  Many 
authors  agree  on  the  main  reason  known  as  the  hindrance  of 


3  Generally  liquid  water  is  not  present  in  bulk  GDL  because  of  its  high 
hydrophobicity. 

4  Bulk  obstruction  of  MPL  is  neglected,  because  its  saturation  is  generally  much 
lower  than  GDL  one  [36].  Anyway  the  GDL  diffusivity  reduction  includes  also  the 
eventual  MPL  one. 
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Table  2 

Parameters  utilized  for  the  calculation. 


Dioues  1.59  X  10-1  X  (T/333)1'75  cm7  cj 

10-5-4163-999.778,rxl04  £ 

DmHffl"oH  C,  X  exp(C2  X  (1/303  -  1/7))  cm2  s"1 

An,H20  4.17  X  10  4  X  22  X  cm2  s_1 

(161exp(— 22)  +  1 )  x  exp(— 2436/7) 
hconv  (@333  K)  0.115  X  102  cm  s_1 

hconv  (@353  K)  0.125  x  102  cm  s  1 

C?ef  1  x  10“3  mol  cm  3 

Cef  7.25  x  10  6  mol  cm  3 

7a  h 

yc  i 

ac  0.6 

of  h 

i*  u  X  exp (75/R  X  (1/353  -  1/T))  A  cm"2 

1.37  x  10  4  x  exp  A  cm  2 

(50,750/R  x  (1/353  -  1/T)) 

0.073  x  exp(1268  x  (1/298  -  1/7))  IT1  cm-1 
Km  W,  m2 


c|F 

ftto 


£o 


Acell 

ha 


6700 

W2  x  (p£0  -p^0) 

W4  x  exp(1029  x  (1/333  -  1/7)) 
C3  x  exp(C4  x  (1/333  -  1/7)) 
(-0.0028  x  (T  -  273)2  -  0.1757  : 
(T- 273) +  1003.8)  x  10-3 
3.56  x  10“4 

2.2  x  exp(5200  x  (1/T  -  1/298)) 
3.5  x  exp(2400  x  (1/T  -  1/298)) 


8314 

0.018 

0.005 

0.025 

0.03 

22.1 

0.08 


Pas 


V 

C  mol1 
J  mor1  K-1 


[43] 

Calibrated 

[64] 

Calibrated 

Calibrated 

[27] 


Calibrated 

[30] 

[3rb  d 


Calibrated 

[38] 


[57] 

Calibrated 

[55] 

Calibrated 

Calibrated 

Calibrated 

[64] 

[64] 

[64] 

[64] 


The  model  is  calibrated  with  respect  to  three  different  typologies  of 
measure:  performance,  water  transport  and  methanol  cross-over. 

The  residuals  between  model  results  and  experimental  data 
have  reasonably  normal  distribution  and  they  are  generally  lower 
than  measurement  uncertainty  for  all  the  three  typologies  of 
measure.  This  implies  that  the  model  is  sufficiently  accurate  to 
estimate  the  available  experimental  results  and  that  the 
phenomena  neglected  with  the  proposed  assumptions  have  minor 
effects.  A  comparison  between  simulated  results  and  experimental 
measurements  of  water  transport  in  two  very  different  operating 
conditions  is  reported  as  example  in  Fig.  7. 

All  the  obtained  values  of  fitting  parameters  are  coherent  with 
those  reported  in  literature.  The  fitting  parameters  of  performance, 
Table  3,  are  composed  of  quantities  affected  by  high  uncertainty, 
demonstrated  by  a  high  variability  in  the  literature  [30,56-58],  The 
effective  methanol  diffusivity  in  the  gas  phase  through  the  diffu¬ 
sion  layer  is  noticeably  uncertain,  considering  the  contribution  of 
CO2  convective  transport,  the  variation  in  void  fraction  and  the 
methanol  mixture  phase  transition.  The  obtained  value  is  very  close 
to  that  reported  in  Ref.  [59],  The  order  of  the  overall  anode  reaction, 
ya,  in  literature  varies  from  0  [56]  to  0.5  [30]  to  1  [60];  the  value 
obtained  by  calibration,  0.065,  evidences  a  very  weak  current 
dependence  on  methanol  concentration.  With  regard  to  the  fitting 
parameters  for  flooding  effects,  the  value  of  Bj  is  very  similar  to  that 
reported  in  Ref.  [43],  instead  for  the  values  of  B2  and  B3  there  is  no 
reference  for  comparison.  Nevertheless  the  reduction  of  cathode 
diffusivity  due  to  liquid  water  flux  through  cathode  GDL  can  be 
considered  reasonable  and  is  discussed  in  the  next  section. 


GM,  1  bar). 


The  fitting  parameters  for  water  transport  are  reported  in 
Table  4.  For  the  fitting  parameter  W2  there  is  no  reference  in 
literature;  anyway  the  resulting  values  of  C™’q,l,  in  the  investigated 
conditions,  are  always  coherent  with  a  highly  hydrated  membrane. 
Electro-osmotic  drag  coefficient  has  a  considerable  variability  [52]: 
most  of  the  works  reports  a  value  of  about  3  at  333  K  [12,19,54], 
some  studies  higher  values,  near  6  [52,63],  moreover  this  value  can 
increase  considerably  when  methanol  is  present  [63],  The  resulting 
value  of  4.6  at  333  K  is  thus  considered  acceptable. 

The  fitting  parameters  reported  in  Table  5  characterize  meth¬ 
anol  cross-over  through  the  membrane.  C3  is  methanol  drag  coef¬ 
ficient,  the  fitting  value,  1.82,  lower  than  water  drag  coefficient, 
confirms  the  influence  of  methanol  concentration  gradient  in  the 
membrane,  as  reported  in  Ref.  [38],  The  dependence  of  electro- 
osmotic  drag  on  temperature,  C4,  is  lower  than  that  reported  in 
Ref.  [37],  evidencing  a  limited  effect  of  temperature. 

Summarizing,  the  values  of  the  parameters  obtained  by  cali¬ 
bration  result  are  acceptable  and  coherent  with  those  reported  in 
literature.  This  consideration,  together  with  the  high  accuracy  of 
the  model  in  reproducing  the  experimental  results,  provides  a  first 
validation  of  the  proposed  interpretation  and  model  of  water 
transport  in  DMFC. 


4.4.  Modeling  results  discussion 

Fig.  8a  illustrates  the  water  fluxes  through  the  GDL  at  different 
current  densities:  when  liquid  permeation  occurs,  the  magnitude 
of  diffusion  flux  diminishes  due  to  the  increased  water  concentra¬ 
tion  in  the  cathode  channel.  At  high  current  densities,  when 
cathode  flooding  occurs,  the  permeation  is  the  predominant  water 
transport  mechanism. 


Table  3 

Fitting  parameters  of  performance. 


Present  work 


7,  cm2  s-1  2.17  X  10-2 

h  -  6.5  x  nr2 

h  -  0.5 

14  A  cm-2  2.69  xio-4 

15  J  mol-1  66,712 

-  37 

B2  -  5.24  x  102 


2.36  x  10-2  [59] 

0  [56],  0.5  [30],  1  [60] 

0.5  [36] 

94.25  x  10-4  [56]  to  2  x  10-4  [36] 
70,000  [61] 

[43] 


Table  4 

Fitting  parameters  of  water  transport. 
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W,  m2  3  X  1(T19  2  10  18  [36]  to  1  x  10~21  [56] 

w2  -  in  x  icr6 

VV3  cm2  s  1  2.38  x  10  2  1.55  x  102  [62] 

VV4  -  4.6  2.9  [54]  to  6  [63] 


Fig.  8b  reports  the  water  fluxes  through  the  membrane.  It  is 
interesting  to  note  that  permeation  flux  is  negligible  and  is  directed 
from  cathode  to  anode.  On  the  contrary  the  diffusion  flux  is  always 
directed  from  anode  to  cathode;  at  low  current  densities  it  is  the 
predominant  water  transport  mechanism,  due  to  a  low  water 
concentration  at  cathode  side,  while  at  high  current  densities  it 
becomes  negligible.  The  resulting  trend  of  membrane  water  fluxes 
is  coherent  with  that  reported  in  Ref.  [36], 

In  this  paragraph  a  specific  analysis  of  flooding  effects  on 
performance  and  a  further  validation  of  the  model  on  MEA  GG 
experimental  results  are  reported. 

4.4.1.  Flooding  effects 

The  developed  model  appears  able  to  predict  the  flooding  onset 
and  to  quantify  its  effects.  Fig.  9  reports  a  comparison  between 
simulation  results  with  and  without  the  implementation  of  flood¬ 
ing  effects,  in  the  same  operating  conditions  of  Fig.  7.  At  lower 
methanol  concentration,  the  water  flow  at  cathode  outlet  is  low 
and  the  modeling  results  with  and  without  the  implementation  of 
flooding  effects  are  superimposed:  this  indicates  that  flooding  does 
not  occur.  Instead  at  higher  methanol  concentration,  where  flood¬ 
ing  is  expected  as  discussed  in  Section  3.1,  a  considerable  difference 
omitting  flooding  model  is  evident  at  high  current  density:  this 
indicates  that  flooding  occurs  significantly  at  current  densities 
higher  than  0.4  A  cm  2. 

A  specific  analysis  is  carried  out  to  evaluate  the  contributions  of 
the  two  flooding  phenomena  considered  in  the  present  model, 
superficial  and  bulk  pore  obstructions.  Fig.  10  illustrates  a  compar¬ 
ison  of  modeling  results  excluding  superficial  or  bulk  pore 
obstruction  correlations,  in  three  different  operating  conditions,  in 
order  to  distinguish  and  quantify  their  effects: 

•  in  Fig.  10a,  the  effect  of  bulk  obstruction  is  much  more  relevant 
than  the  superficial  one; 

•  in  Fig.  10b,  water  concentration  in  cathode  channel  increases 
due  to  the  lower  airflow  rate,  enhancing  water  condensation 
and  superficial  obstruction,  nevertheless  bulk  obstruction  is 
still  more  influent; 

•  in  Fig.  10c,  the  saturation  concentration  diminishes  due  to  the 
lower  temperature,  leading  to  a  further  enhancement  of  water 
condensation  on  GDL  surface.  In  this  case  bulk  and  superficial 
obstructions  are  comparable. 

Generally  the  effect  of  bulk  obstruction  is  more  relevant 
compared  to  the  superficial  obstruction  one.  However  both  types  of 


Table  5 

Fitting  parameters  of  methanol  crossover. 

Present  work  Literature 

~  cm2  s-1  2.17  x  to-6  0.5  x  1CT6  [23]  to  4.9  x  1CT6  [36] 

C2  K  1976  2436  [36] 

C3  -  1.81  2.9  [37]  to  0.87  [38] 

C4  K  410  1029  [37] 


obstruction  must  be  taken  into  account  to  accurately  reproduce 
experimental  data:  calibrating  the  model  without  one  type  of 
obstruction,  the  model  accuracy  decreases  dramatically. 

Fig.  11  reports  the  current  density  profiles  along  the  channel  at 
different  voltages.  From  Fig.  9  it  is  possible  to  figure  out  that 


CM,  1  bar). 
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flooding  occurs  at  voltages  lower  than  0.2  V.  Analyzing  the  current 
density  profiles,  a  flooding  effect  is  already  evident  at  0.2  V:  the 
current  density  profile  presents  a  more  steep  decrease  than  at  0.3  V. 
The  slope  decreases  even  more  at  0.1  V,  evidencing  a  hindered  fuel 
cell  operation  near  channel  outlet,  caused  by  low  oxygen  concen¬ 
tration  at  cathode  electrode.  Such  high  current  density  gradient  not 
only  worsens  performance,  but  could  also  determine  inhomoge¬ 
neous  and  aggravated  components  degradation. 

Fig.  12  illustrates  the  reduction  of  cathode  GDL  effective  diffu- 
sivity5  at  different  voltages.  At  0.3  V  the  reduction  is  in  the  order  of 
20%  and  it  is  only  due  to  liquid  permeation,  i.e.  bulk  pore 
obstruction.  At  0.2  V  the  reduction  is  more  pronounced,  in  agree¬ 
ment  with  the  higher  water  production  and  the  steeper  current 
density  profile,  Fig.  11.  In  this  operating  condition  the  water 
concentration  in  the  channel  reaches  the  saturation  value  approx¬ 
imately  at  two-thirds  of  the  channel  and  consequently  the  effective 
diffusivity  presents  a  further  reduction  due  to  superficial  pores 
obstruction.  At  0.1  V  the  superficial  pores  obstruction  occurs  nearly 
in  the  middle  of  the  channel  and  the  overall  diffusivity  at  channel 
end  is  reduced  by  one  order  of  magnitude:  in  this  condition  the 
cathode  GDL  is  severely  flooded. 

4.4.2.  Cathode  MPL  influence 

The  model  is  further  validated  on  the  experimental  data  of  the 
MEA  GG  (66  measurement  points,  three  types  of  measures),  Table  1. 
Two  parameters,  affected  by  the  absence  of  cathode  MPL,  are  cali¬ 
brated.  The  parameters  reported  in  Tables  2—5  remain  valid,  with 
the  exception  of  the  five  parameters,  related  to  cathode  GDL, 
reported  in  Table  6. 

It  is  interesting  to  note  that  the  fitting  parameter  B2,  that 
quantify  the  bulk  obstruction  effect,  is  higher  than  the  MEA  GM 
one,  Table  3:  the  bulk  obstruction  effect  appears  more  relevant 
without  cathode  MPL  at  constant  water  permeation. 

The  fitting  parameter  W2,  that  characterizes  the  dependence  of 
C™2’o’l  on  capillary  pressure,  is  higher  compared  to  MEA  GM  one, 
Tafele  4:  the  lower  GDL  hydrophobicity  implies  a  lower  capillary 
pressure  at  constant  membrane  water  concentration  [17],  The  pum 
value  is  assumed  lower  according  to  a  less  hydrophobic  GDL 
compared  to  MEA  GM  one  [55]. 

Although  the  uncertainty  of  material  characteristics  and  the 
unpredictable  interactions  among  components,  the  model  is  able  to 
reproduce  accurately  the  experimental  data  of  MEA  GG,  Fig.  13. 
Moreover  the  model  residuals  are  close  to  measurement  uncer¬ 
tainty  and  the  fitting  parameters  assume  reasonable  values.  This 
result  provides  a  further  validation  of  both  water  transport  and 
flooding  proposed  models. 

The  analysis  of  modeling  results  permits  also  to  evaluate  the 
cathode  MPL  influence.  Fig.  13  points  out  the  effect  of  flooding  in 
two  considerably  different  operating  conditions.  The  behavior  is 


5  The  GDL  +  MPL  effective  diffusivity  is  calculated  with  Formula  (16)  [43]. 


Fig.  11.  Current  density  profiles  at  different  voltages  (MEA  GM,  met  6.5%,  353  K,  air 
1.14  g  min-1, 1  bar). 


similar  to  that  reported  in  Fig.  9.  At  higher  methanol  concentration, 
two  differences  are  appreciable:  flooding  effect  magnitude 
increases  and  flooding  onset  occurs  at  lower  current  densities, 
nearly  at  0.35  A  cm-2.  These  aspects  are  coherent  respectively  with 
the  increase  of  bulk  obstruction  effect,  B2,  and  the  reduction  of 
breakthrough  pressure,  pum-  This  analysis  confirms  the  explanation 
suggested  discussing  Fig.  4:  the  presence  of  cathode  MPL  increases 
mass  transport  resistance,  reduces  considerably  water  loss  toward 
cathode  channel  and  increases  DMFC  performance  exclusively  in 
case  of  significant  GDL  flooding. 


Channel  lenght  [%] 

Fig.  12.  Cathode  GDL  +  MPL  effective  diffusivity  profiles  at  different  voltages  (MEA 
GM,  met  6.5%,  353  K,  air  1.14  g  min1, 1  bar). 
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(MEA  GG,  1  bar). 


Table  6 

Fitted  and  assumed  parameters  for  model  validation  (MEA  GG). 

B2  -  6.21  x  102 

W2  -  2.76  x  10  5 

(mpl  cm  0 

/gdl  cm  0.03 

Plim,GDL  Pa  1700  [55] 


5.  Conclusions 

Water  transport  and  flooding  in  DMFC  are  investigated  through 
both  experimental  and  modeling  analyses.  An  exhaustive  inter¬ 
pretation  of  the  mechanisms  that  regulate  these  phenomena  is 
proposed.  The  implemented  model  is  able  to  reproduce  accurately 
the  experimental  data.  Such  model  is  validated  by  means  of  three 
sets  of  experimental  data  (polarization  curve,  methanol  crossover, 
water  flow  at  cathode  outlet),  in  an  extensive  range  of  operating 
conditions.  The  model  includes  also  two  correlations  to  reproduce 
the  effects  of  cathode  GDL  flooding.  A  further  validation  of  the 
model  is  carried  out  on  the  experimental  data  regarding  a  second 
DMFC  with  a  different  cathode  GDL 

The  modeling  analysis  permits  also  to  examine  both  onset  and 
magnitude  of  cathode  GDL  flooding  and  the  influence  of  the  MPL 
presence.  The  main  conclusions  on  water  transport  and  flooding 
phenomena  are  the  followings: 

•  water  transport  through  cathode  GDL  is  regulated  by  two 
mechanisms:  vapor  diffusion  and  liquid  water  permeation; 
analyzing  water  flow  at  cathode  outlet,  the  first  determines 
a  plateau  at  low  current  density;  the  second  a  linear  trend  at 
high  current  density;  influence  of  temperature,  pressure  and 
airflow  confirm  the  expectations; 

•  water  diffusion  through  cathode  GDL  is  regulated  by  vapor 
concentration  gradient;  in  most  of  the  investigated  conditions 
vapor  concentration  in  the  electrode  is  close  to  saturation 
value,  in  agreement  with  a  highly  hydrated  membrane; 

•  liquid  water  permeation  through  cathode  GDL  occurs  when 
water  pressure  exceeds  a  threshold  value,  related  to  GDL 
characteristics;  the  permeation  linear  trend  with  current 
density  is  due  to  electro-osmotic  drag  and  water  production  at 
cathode  electrode; 

•  liquid  water  permeation  through  the  membrane  is  directed 
from  cathode  to  anode,  while  liquid  water  diffusion  has  the 
opposite  direction  with  generally  higher  values;  at  low  current 


densities  diffusion  is  the  predominant  water  transport 
mechanism; 

•  to  reproduce  the  effects  of  cathode  GDL  flooding,  two  mecha¬ 
nisms  have  to  be  considered  simultaneously,  giving  compa¬ 
rable  contributions:  superficial  and  bulk  pore  obstructions;  the 
first  is  proportional  to  liquid  water  concentration  in  cathode 
channel;  the  latter  is  proportional  to  liquid  water  permeation; 

•  a  correlation  to  reproduce  bulk  pore  obstruction  is  proposed 
for  two  GDLs,  with  and  without  MPL;  the  magnitude  of  bulk 
obstruction  effect  is  more  relevant  without  cathode  MPL; 

•  the  addition  of  the  MPL  on  cathode  GDL  increases  the  mass 
transfer  resistance,  causing  a  general  reduction  of  water 
transport  through  cathode  GDL;  the  performance  decreases 
due  to  lower  oxygen  concentration  when  no  flooding  occurs, 
while  performance  increases  with  severe  flooding,  because  the 
MPL  limits  its  effect. 
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List  of  symbols 


C  species  concentration  in  channel  (mol  cm'3) 

C  time-average  concentration  in  channel  (mol  cm-3) 

Cref  reference  concentration  (mol  cm-3) 

Deff  effective  diffusivity  in  diffusion  layer  (cm2  s”1) 

Dm  effective  diffusivity  in  membrane  (cm2  s'1) 
bconv  convective  mass  transport  coefficient  (cm  s'1) 

Eo  ideal  potential  difference  (V) 

F  Faraday  constant  (C  mol'1) 

h  channel  height  and  width  (cm) 

i  local  current  density  (A  cm'2) 

f  exchange  current  density  (A  cm'2) 

k  Tafel  constant  (A  cm'2) 

f<H  Henry  constant  (mol  J) 

L  channel  length  (cm) 

lm  membrane  thickness  (cm) 

1  diffusion  layer  thickness  (cm) 

M  molecular  weight  (g  mol'1) 

m  inlet  mass  flow  rate  (g  min'1) 

Km  membrane  permeability  (m2) 

JVH2°  water  flux  (mol  cm'2  s'1) 

X  molar  or  mass  fraction  (%  or  wt%) 

N  flux  (mol  cm”2  s”1) 

rid  water  drag  coefficient 

n dx  methanol  drag  coefficient 

p  channel  pressure  (Pa) 

Psat  saturation  pressure  (MPa) 

Piim  threshold  breakthrough  pressure  (Pa) 

R  universal  gas  constant  (J  mol'1  K”1) 

Rc  mass  transport  coefficient  in  gas-phase  (s  cm”1) 

T  fuel  cell  temperature  (K) 

v  local  velocity  in  channel  (cm  s'1) 

x  x  coordinate  (cm) 

Creek  symbol 

a  Tafel  transport  coefficient 

y  reaction  order 

e  volumetric  void  fraction 

t]  polarization  (V) 

<rm  membrane  conductibility  (O'1  cm”1) 

p  compound  density  (g  cm”3) 
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/i  compound  viscosity  (Pa  s) 

Superscript 

a  relative  to  anode 

c  relative  to  cathode 

t  relative  to  catalyst  layer 

G  relative  to  gas  phase 

L  relative  to  liquid  phase 

sat  relative  to  saturation 

Subscript 

GDL  relative  to  gas  diffusion  layer 

MPL  relative  to  micro  porous  layer 

m  relative  to  membrane 

CO2  relative  to  carbon  dioxide 

H2O  relative  to  water 

CH3OH  relative  to  methanol 

N2  relative  to  nitrogen 

O2  relative  to  oxygen 

air  relative  to  air 

in  relative  to  channel  inlet 

out  relative  to  channel  outlet 

drag  relative  to  electro-osmotic  drag 

diff  relative  to  diffusion 

cross  relative  to  cross-over 

perm  relative  to  permeation 

ox  relative  to  oxidation 

rd  relative  to  reduction 

tot  relative  to  total 

conv  relative  to  convective 
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